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Abstract 



In this work we present a model for the determination of the interaction energy for triplet and singlet 
states in atoms with incomplete filled shells. Our model includes the modification of the Coulomb's 
law by the interaction between the magnetic moments of the electrons and a Heisenberg term. We 
find that the energy of the triplet state is lower than the energy of the singlet state. We calculate the 
interaction energy between the electrons from the adjacent atoms in fee and bec lattices and we find 
that the minimum interaction energy is attained for the triplet state. The result is presented for the 
interaction between the electrons of the first coordination group. The interaction energy which aligns 
the spins is used to evaluate the Curie temperature in a mean field model. Compared to previous 



• models, our simple model fits the experimental data. 
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£cj ■ 1 Introduction 
in 

Various models and approximations are used to explain the emergence of the ferromanetic state. The 
physical effect that produces magnetic ordering of adjacent magnetic moments is the same that leads 
to magnetic ordering within an ion and produces Hund's rules. Interactions between spins are actually 
— , ' just a convenient way to record the end result of electrostatic repulsion pQ. Two electrons circling a 
nucleus reduce their Coulomb interaction by adopting an antisymmetric wave function that vanishes 
whenever they come near each other. The Pauli principle demands that the overall wave function must 
be antisymmetric, so the spin wave function must be symmetric. The electrons lower their energy by 
adopting the same spin state and developing a local magnetic moment. The important concepts were 
introduced by Heisenbergj2] and Stoner[3]. It was shown by Hubbard[4j that the correlation effects 
will lower the energy of non-magnetic states more than that of the ferromagnetic states, and so make 
the condition for ferromagnetism more stringent. The spin wave excitations give a way to explain the 
temperature dependence of the magnetization in magnetic materials at low temperatures [5j. 

In this paper, using an effective field Hamiltonian (detailed in Ref.[6]), we investigate the influence 
of the magnetic moments on the electron-electron interaction. We first discuss the Heisenberg term 
that emerges in the energy of interaction. Next, we present the modification of Coulomb's law by the 
interaction between the magnetic moments of the electrons. Further, we discuss the case of an isolated 
atom with applications to Hund's rules and the condition for ferromagnetism in fee and bec lattices. 
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2 Model and formalism 



In this section we present the electron-electron interaction by using an effective Hamiltonian[6] in which 
we separate explicitly the Heisenberg term. In the interaction picture the effective Hamiltonian of the 
electron-electron interaction is given by the expression (see Appendix A): 
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where I? is a coupling constant, m is the mass of an electron, r is the distance between the two electrons, 
po is the massive density of the interacting field, Dr/cP is the "massless density" of the interacting field, 
Uq=cq is the classical oscillation frequency of the interacting field, w q o is the oscillation frequency of 
an electron, q is the wave vector of the interacting field, qo is the wave vector of the boson associated 
with the electron, k is the wave vector of the electron and = hk 2 /2m. We introduce the creation and 
annihilation operators Cj a + and Cj a of an electron of spin a in the orbital state /( r -Tj): 
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and the spin operators: 



where a a (a = x, y, z) are the Pauli matrices. 
After some algebra we obtain: 
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with nj = c ^j U c 3(T the occupation number operator. 
Therefore we can rewrite Eq.([T]) as: 
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The first term in Eq.@ is a "pure" Coulomb term while the second term represents the Heisenberg term: 

H H = 2Y^K jf s j -s j/ (6) 

33' 

3 Effect of the magnetic moments on the interaction energy 

We first consider two electrons situated at a distance r = |r, - rj/\. The expectation value of the first 
term from Eq.(j3]) is: 

E c = \^K jj ,n j n r (7) 
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For n q = n q Q = and rij = ny = 1, we have: 
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where in the presence of the magnetic field we introduce the vector potential and thus we substitute qo- 
r by qo- r - (e/h)§ A • dr. 

The energy levels of atomic electrons are affected by the interaction between the electron spin magnetic 
moment and the orbital magnetic moment of the electron (spin-orbit coupling). It can be visualized as a 
magnetic field caused by the electron's orbital motion interacting with the spin magnetic moment. The 
effective magnetic field can be expressed in terms of the electron orbital angular momentum. We consider 
the vector potential A = ^^r" where fi is the magnetic dipole moment and r is a vector from the middle 
of the loop to an observation point. An electron in a stationary state in an atom, having a definite 
angular momentum projection L z = hmi (mi the quantum magnetic number), has a magnetic moment 

= HB m i where ^b = eh/2m is the Bohr magneton. The theory and experiments demonstrate that 
the free electron has a magnetic moment equal to the Bohr magneton [ib , and a spin momentum s, the 
projection of which on a specified direction are s z = ±h/2 = hm s where m s = ±1/2 is the spin quantum 

(s) 

number. For fj,\ = ^B9 rn s with g = 2 we obtain: 
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where (h/e) ran and (h/e) m s j are the flux vectors. For q$j = q$ji =q$ , we have: 
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In these conditions q = q — ^- ( m '+ 2ms ) x where x is a unit vector perpendicular to r and //. 
When the interacting field is a photon field, then ro = 0. For a quasi free electron, — £k-q << <^qj 



as: 



%q^j2m (with m an effective electron mass). 
To calculate (Eq.(|5|)) we consider the two sums separately. The first sum can be easily calculated 
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where we have considered Q = 4irr 3 /3. The second sum can be calculated as: 
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In these condition Ec (Eq.(ff|)) becomes: 
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The upper limit of go as 0.767r/r is imposed by the constraint condition 
m s j = 1/2, Tq reduces to: 
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where ro =rj = ry is the radius of the electron orbit and r = rjji. For mij — m\y = 1 we get: 
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For m s j = 1/2, m SJ ' = -1/2, we obtain: 
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For mij = mij' the above relation reduces to: 
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We now calculate the second term of the interaction energy (Eq.Q). For an interacting pair of spins 
s\ = s 2 = 1/2, the expectation value of the Heisenberg term (Eq.flSJ)) is written as: 
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with S = si +s 2 . In the singlet state, S = and Eh= - (3/2) Ec- In the triplet state, S 
(l/2)Ec, where Ec is given by Eq.([9|). The spin dependent energy of interaction is: 
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The term responsible for magnetic ordering is the second term: 
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For the singlet state (5 = 0), we get Ef,j= 0, while for the triplet state (5 = 1) we obtain: 
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We find that for cos To < 0, we have E l M < E S M and the interaction energy for the triplet state is 



lower than that associated with the singlet state. Moreover, for an average value of r equal to r Q , Y s ^ n9let 
_ 2p^ip et j n FigJJJ we /\j a — cosTq — cos(2rg), that is the difference between the function for 
the triplet state and the function for the singlet state. It is observed that for To = vr the amplitude of 
this function has an absolute maximum, which is negative and therefore the interaction energy Eq.(jl]) 
for the triplet state is lower than that associated with the singlet state. Thus, for an atom with two 
valence electrons in a shell with I > 1, the interaction energy between the two electrons when they occupy 
different L z states (all degenerate in energy) and adopt a triplet spin state, is smaller than the interaction 
energy when they occupy the same L z state, in which case they must have opposite spins. Therefore, 
the triplet state is favorable. This could explain Hund's rules related to the organization of electrons in 
incompletely filled shells. Finally, we note that if we do not separate the Heisenberg term, and do not 
take into account the magnetic moments, the expectation value of Hj (Eq.(JT])) would be[6 : 
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which is just the Coulomb's law where the constant 1/147 is close to the fine structure constant, a 
1/137. This difference appears due to constraint in the evaluation of the integral over go- 



4 Results for fee and bee lattices 

We analyze now the case of the face-centered cubic (fee) and base-centered cubic (bec) lattices and we 
determine the conditions for the emergence of ferromagnetism. First, we consider that at some instant 
a given atom has its total spin in the up direction. The intra-atomic interactions are of such a nature 
that this atom tends to attract electrons with spin up and repel those with spin down. In the case of 
the fee lattice every atom is surrounded by 12 nearest neighbors, which have coordinates (a/2)(±l,±l,0), 
(a/2)(±l,0,±l), (a/2)(0,±l,±l). We consider that the magnetic moment is /^(0,0,//). Then, the vector 
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Figure 1: The difference between the functions for the triplet and singlet states A/ , as a function of G Q , 
in an isolated atom with an incompletely filled shell. 



potential can be chosen in the following form A(A x ,A y ,0). For rij 
the sum in the interaction energy (Eq.([T])) becomes: 
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where // cc (ro) = (cos To + sin To) 2 + 2 (cos To + sin To) and To is given by Eq. pT]) for the triplet state 
and by Eq. (|13p for the singlet state. The Coulomb interaction is then given by: 
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The term responsible for magnetic ordering is: 
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which has a minimum for To = 2nir — tt/4. 

In a bec structure, every atom is surrounded by 8 nearest neighbors which have the coordinates (a/2)( 
± 1, ± 1, ± 1 ). the sum in the interaction energy (Eq.([T])) now becomes: 
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with /fe cc (r ) = (cosTo + sinTo) 2 . 

The energy of interaction is given by the expression: 

he 



Ei 



1 q (cos r o + sin r ) 2 

IT 6 



18tt r 

and the term responsible for magnetic ordering is: 



5 + + 



4he 

e m = -^fbccS (S + l) = -AEjS (S + 1) (21) 
which is somehow different from the energy Eq. (|'20p obtained for the fee structure. 

4.1 Discussion 



From Eq. (|20p and f)21 1) it appears that there is an oscillating interaction between the two electrons belong- 
ing to different atoms. The nature of this interaction may be different than that which is encountered in 
the case of interaction between two ions by indirect exchange [7J. In the later case, the electrons belonging 
to an ion flip a conduction electron which then travels to another site and interacts with the spin of the 
ion of the second site (RKKY interaction). Perhaps the most significant application of the RKKY theory 
has been to the theory of giant magneto resistance (GMR)|8, 9j. GMR was discovered when the coupling 
between thin layers of magnetic materials separated by a non-magnetic spacer material was found to 
oscillate between ferromagnetic and anti ferromagnetic as a function of the distance between the layers. 
The period of an oscillation is determined by the Fermi wave vector, in the case of free electron gas, via 
A = 27t/2/cf = Tr/fep. The phase of the cosine in RKKY theory is kpr while in our model this phase is 
Tq = ne 2 /mr = l/k^r where k$ = m/ne 2 . In this case, the "period" of an oscillation is dependent of 
r. From the condition — Tq^ 1 ^ = 2ir, results A n = r n — r n+ \ = 2kor n r n+ \. The effective mass of 
the electron in the rf-band is larger than the free electron mass. For r n = 3. 3 A and m = 5m-o results 
A n ^ OAA a value which is larger than the Compton wavelength, h/mc = 8 x 10~ 3 A, and is lower than 
the period of the RKKY oscillation, which is of the order of 3A. The oscillatory behavior is the result of 
the interference of two oscillating fields generated by the magnetic moments of the two electrons. Each 
fringe arises from a definite difference in phase. 

The Curie temperature may be estimated as in the mean field theory applied to the spin-wave model[10, 

ITS: 

where M is the magnetic moment per atom, is the Bohr magneton, _E q is the spin- wave energy which 
is related to the exchange parameter Jqj by a simple Fourier transformation: 



E ^ = ~w £ Joj ^ ~ exp ' R ^ 



Substituting the exchange parameter by our interaction parameter used in Eqs. (l20p - (|2ip finally we 
obtain: 

k B T c = ^AEj (22) 



where AEj is given by the above expressions. For Co, Tc Xp = 1390 K, we obtain from Eq. (120p and Eq. (|22j) . 
ffee = 0.064 which is appropriate for an effective mass of the electron in (f-band m = 5.67mo- For Ni, 
Tc xp = 630 K, we obtain f fcc = 0.029 and m = 5.39m . For Fe, T e c xv = 1045 K, we obtain f bcc = 0.037 
and m = 5.71mo- 

A more elaborate method is the random phase approximation [10 1 E l H^l ■ Also, the Curie temperature 
can be estimated numerically by employing the method of Monte Carlo simulations applied to effective 
Heisenberg Hamiltonian|13j . 
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5 Conclusions 



Using an effective field Hamiltonian (Ref.[6]) we have shown that the Coulomb's law is modified by 
the spin-spin interaction. First, we separate the Heisenberg term in the energy of interaction. Next, 
we take into account the effect of the spin magnetic moment. The Coulomb interaction is modulated 
by a cosine term whose argument depends on the spin magnetic moment and on the distance between 
the two electrons. Evidently, this situation occurs in the absence of the magnetic field, for example in 
superconductors. Further, we have studied the influence of the electron magnetic moments, both orbital 
and spin, on the electron-electron interaction in an isolated atom. We have found that in an incomplete 
shell the interaction energy of the triplet state is smaller than the energy of the singlet state. On this 
basis may be explained Hund's rules. The condition of ferromagnetism is studied in both fee and bee 
lattices of transition metals. By using the mean field approximation we have estimated the effective mass 
of the electron in e?-band, by knowing the Curie temperature. The coupling energy, due to the magnetic 
field generated by the two interacting electrons, oscillate with a "period" l/kor where ko = va/ire 2 . This 
period of oscillation is dependent of r and for r = 3.3A and m/rriQ =5, is equal to 0.4 A, a value which 
is smaller than the RKKY period which is of the order of 3 A. 



Appendix A 

We provide here a detailed view of the interaction energy as expressed by Eq.([T|). 

The electron will always carry with it a lattice polarization field. The composite particle, electron 
plus phonon field, is called a polaron. It has a larger effective mass than the electron in the unperturbed 
lattice. By analogy, in our model, we consider a coupling between an electron and a boson. 

Consider a linear chain of N bodies, separated by a distance R. The Hamiltonian operator of the 
interacting bodies (electrons) and the boson connecting field takes the general form: 



H — H 0fi i + H 0iP h + Hj 

where 

k,cr 

is the Hamiltonian of the electrons of mass m, , Cfc i(J are the electron creation and annihilation operators, 
k is the wave vector of an electron and s is the spin quantum number, 
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where w q = \ ° t+D J' q J \ s the classical oscillation frequency, a is the restoring force constant, D is the 

coupling constant, p is the linear density of flux lines, the boson creation and annihilation 

operators, p = po + Dr/c 2 , where po is the density of the interacting field, if this is a massive field, c is 
the velocity of the boson waves. The interaction Hamiltonian operator Hj is given by the expression: 
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6q , 6q are creation and annihilation operators associated with the electron oscillations, e q i denotes the 
polarization vectors and xi a ) is the s pi n wave function. s n is the displacement of a body near its 
equilibrium position in the direction of r and, in the approximation of nearest neighbors, it is assumed 
that D does not depend on n. The Hamiltonian of interaction between electrons via bosons becomes: 

1 D fs—^ ^ H f u + \ frq • q' 

1 ~ ~W2^ J ^ ^ 2^Z~ V &q « + b -*>) ~2^Z 

J n k,k',q,q',q ,a,a' q ° ' Q 
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q, q' are the wave vectors associated with the bosons of the connecting field, q Q is the wave vector 
associated with the oscillations of the electron, and k, k' are the wave vectors of the electrons. Consider 
the integral over z: 

J g-ife+g'+k-k'+qo)-*^ = NrA ( q + q / + k _ k / + qo ) ( 23 ) 

where A(x) = 1 for x = and A(x) = 0, otherwise. In the crystal bulk Nr is replaced by V = where 
Q is the volume of a unit cell and N is the number of unit cells. We obtain: 
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Approximating 

(a q + atq) (a-q' + « q /) ~ a q a q / + al q «- q ' 
and choosing in Eq. (|23p . q' =qo and k' = k + q, the effective Hamiltonian is given by: 
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The expectation value of the energy of Hj{ is the energy of the electron-electron interaction which is 
presented in the text. The expectation value of the energy of U e JJ is the self energy of the electron and 
is used to calculate, for example, the polar on energy[B]. 
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